Estimate of the weight of empty space based on astronomical 
observations. 
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Abstract - As a consequence of the equivalence principle and of the existence of a negative 
vacuum energy, we show how a weak but universal linear response of the vacuum to a local 
gravitational potential suffices to explain the main features of so-called Pioneer anomalies as well 
as the apparent departure from the third Kepler law which has been observed at the level of 
numerous galaxies. 



Introduction. — It is nowadays widely accepted that 
it is impossible to explain the deviations from Kepler's 
third law exhibited by stars at the level of numerous galax- 
ies without formulating extra- hypotheses. The two main 
hypotheses that are presently [Tl[2] competing are (1) the 
MOND theory and (2) the dark matter hypothesis. We 
propose a third approach (VIND for Vacuum Influenced 
Newtonian Dynamics) in which Newtonian dynamics is 
preserved while, in the absence of dark matter, it is the 
linear response of the vacuum energy which is supposedly 
responsible for apparent deviations to Kepler's third law. 
The paper is structured as follows. In the first section, we 
introduce the MOND, dark matter and (3) VIND mod- 
els. In the second section we investigate the modifications 
to standard Newtonian dynamics deriving from our VIND 
model at small distances compared to the typical scale of 
the model (L). This allows us to estimate the value of L, 
on the basis of so-called Pioneer anomalies. We then in- 
vestigate the modifications at large distances and obtain a 
good agreement between our predictions and MOND phe- 
nomenology. The last section is devoted to open questions, 
discussion and conclusion. 

MOND, dark matter and VIND models. (1) 

In the MOND theory [3] ("MOdificd Newton Dynamics") 
one assumes that Newton's law must be modified in a ad 
hoc fashion whenever we consider very small accelerations 
(relatively yo the inertial frame attached to "fixed" stars). 
It is motivated by the observed violations of Kepler's third 



law, which establishes that, as is the case in the solar sys- 
tem, celestial objects follow trajectories for which "the 
square of the period T varies like the cube of the mean 
distance R to the sun" . In the case of circular orbits, this 
reflects Hookc's law: ^ R = a^*^™*"^ is proportional to 
up to constants. In other words, the Newtonian force 
m • a undergone by a test-particle of mass m at distance 
i? of a ponctual source varies like On the contrary, 

MOND theory in its simplest form predicts that the (in- 
tensity of the) acceleration undergone by a particle test in 
a very weak gravitational field (for instance far away from 
other gravitational sources) obeys 
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where ao is a new universal constant introduced by Mil- 
grom. It is easy to check that, far away from a ponctual 
source of mass A/, MOND theory predicts that 
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-^GnM ■ ao/r^ = -{^JGnM ■ ao)/r. 



(2) 

Therefore, the acceleration undergone by a test-particle 
at distance i? of a ponctual source varies like 1/i? so that 
(^)^i? is attractive and proportional to In other 
words, the velocity R/T is constant (does not depend on 
R) for any circular orbit, which is a fundamental prop- 
erty of MOND theories, that has been sometimes infirmcd 
in the practice but has nevertheles been confirmed by nu- 
merous observations [4l[5] (this is sometimes called the flat 
velocity distribution) . 
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(2) in the dark matter approach it is assumed that the 
departures from Kepler's third law are due to the pres- 
ence of a huge concentration of dark matter, invisible to 
our telescopes and accelerators but responsible for the ap- 
pearance of large accelerations which cannot be explained 
otherwise, if we stick to the Newtonian paradigm. 

Both approaches present difficulties. The MOND ap- 
proach can be considered as a clever and economic phe- 
nomcnological description of what happens in a large num- 
ber of galaxies, but it suffers from a clear theoretical back- 
ground and motivation. For instance it is very difficult to 
conciliate MOND theories and general relativity, the lat- 
ter predicting, roughly, that at slow velocities, Newtonian 
gravitation prevails. It also suffers from clear experimen- 
tal confirmations at the scale of clusters of galaxies and at 
the scale of the solar system. 

The introduction of dark matter can also be considered 
to be an ad hoc hypothesis being given that, so far, all at- 
tempts to find a trace left by this hypothetical dark mat- 
ter outside from the specific context of galactic rotations 
curves failed. At CERN for instance, none exotic particle 
described by dark matter theorists has been observed so 
far. 

(3) The goal of the present paper is to show that the 
MOND predictions (or at least an "average" version of 
them) can be infered from an assumption that has the ad- 
vantage to respect Newton dynamics. The price to pay 
is that we must assume that, everywhere in the cosmos, 
empty space exerts a weak back-reaction onto the gravi- 
tational potential. Otherwise, our results are derived in 
the framework of Newton's theory of gravitation. In par- 
ticular we do not need to modify Newton's dynamical law 
{F = ma) at small accelerations, and the gravitational 
potential still evolves according to Poisson equation 



-2px) + i47rGNC ^p' 



— 2 ^matter 
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A(p = AttGnp, 



(3) 



where is the gravitational potential, Gjv the Newton 
gravitational constant and p the density of mass-energy. 
It has been recognised since several years that in the New- 
tonian limit the density of mass-energy ought to contain 
a contribution due to the so-called cosmological constant, 
which leads to the modified Poisson equation [6] 



A(/) = 47rGjv(p" 



- 2pa), 



(4) 



where px = A/SttGat and A is Einstein's cosmological con- 
stant. It is not completely clear whether this contribution 
is due to quantum vacuum fluctuations or not but here 
we shall take for granted that empty space possesses a 
negative mass density, constant in space. 

The novelty of our approach is that we also assume that, 
as a consequence of the mass-energy equivalence principle, 
the matter and the vacuum do not only contribute to the 
density of mass-energy as isolated source terms, but that 
they also exert a back reaction proportional to the gravi- 
tational field: 



where = SttGnc'^px = Ac^^. 

We are not interested here in computing effects 
due to the linear in back reaction of matter 
(47rGArc~^p™"**'^''0). They can be evaluated indepen- 
dently and in general they only result in a small correction 
to usual Newtonian dynamics. For instance, at the level 
of the solar systems where the gravitational potential is 
everywhere small relatively to c^, they does not modify 
much the Newtonian, zero-order dynamics. Nevertheless, 
the back reaction due to the vacuum energy, despite of 
the fact that it is weak, is present everywhere in space 
and its cumulative effects may not be neglected. Our goal 
is to show how these effects provide a plausible explana- 
tion of the observed deviations to Kepler's third law, and 
of so-called Pioneer anomalies as well... 

In summary. Vacuum Influenced Newton Dynamics de- 
rives from the assumption that Newton dynamics is not 
modified but "only" indirectly infiucnccd by the appear- 
ance of a new source term k'^(p that we ultimately attribute 
to vacuurr0. 

VIND in the near and far field regimes . For- 
mally, the equation , in vacuum, can be solved once we 
know the associated Green's function, which is a solution 
of Green's equation 



AG(x, x') + /c^G(x, x') = (5(x - x'). 



(6) 



where (5(x — x') represents the three-dimensional delta 
function. As is well-known, even if we impose that the 
Green function is isotropic (that is to say, it depends only 
on the distance between x and x'), it is not unique in the 
sense that all Green functions differ by a solution of the 
homogeneous Helmholtz equation. We are interested in 
a real and isotropic Green function, which suggests the 
choice 

— cos(fc(|x — x'l) 



G(x - x') = 



47r|x - x' 
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We are free to add to the Green function defined by 
equation Q a real multiple of the harmonic function 
sm(k(\x~x I) ^ However, in realistic situations, the weight 

47r|x— X I ' ' ^ 

of such a solution remains very small, and the essential 
conclusions of all our analysis remain unchanged. 

In particular, convoluting the omnipresent vacuum en- 
ergy —2/9 A everywhere in space with the Green function 
yields a finite constant additive correction to the poten- 
tial. The contribution in a location x essentially comes 
from a region situated at a finite number of wavelenghts 
(denoted L, L = 27r/fc) around x. It is not necessary to 
make particular assumptions about the structure of the 
universe (e.g. is it finite, open or else) in order to renor- 
malise properly the gravitational potential in the present 
case, excepted that L supposedly differs from 0. We shall 
thus consistently forget in what follows the presence of the 



A(/) = 47rGjv(/5" 



2pa)(1 + 0/c2) 



^See ref. [7| for a similar angle of approach. 



p-2 



Estimate of the weight of empty space 



source term —2p\ in equation ([S]) and focus on the changes 
brought by the hnear vacuum response k'^cj). 

Near field regime. At smaU distances (|x — x'| << L), 
the Green function reduces at the zeroth order to the usual 
Green function of Newton's theory 4^ | ~lx'| ' which means 
that our model is nearly equivalent to Newton's model for 
gravitation at small distances from the sources of matter. 
As it is granted that, in the solar system, Newton's model 
works very welfl we must impose that L is quite larger 
than the size of the solar system. It is however interest- 
ing to take account of slight but systematic deviations of 
Newton's theory which have been observed with slow ob- 
jects (the Pioneer 10 and 11 probes), and constitute what 
is called the "Pioneer anomaly" . 

In short, it has been repeatedly observed at the level of 
these probes (twice) that a constant centripetal accelera- 
tion seemingly prevails in our solar system. This acceler- 
ation is constant in time and space and directed towards 
the sun. Its mag nitude is [8] of the order of 8 • 10-^° m/s^. 

Our model precisely predicts such an effect, as can be 
seen by developing the Green function ^ to the lowest 
order in k which reads 



and equal to < p" 



>sun, it is relatively easy to esti- 



G(x - x') 
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Sufficiently far from the sun (where x — x' is in good ap- 
proximation constant in magnitude and direction for all 
source terms inside the sun (such that |x'| < Rsun)), so 
that we expect that a constant radial and centripetal force 
of intensity Gn ■ Msun ■ k'^/2 must exist, which adds to the 
usual Newtonian solar attraction. Identifying this force 
with the Pioneer anomaly, we derive an estimate of k: 
k^ = 2apioneer/ (Gn ' Msmi)- Injecting in this relation the 

typical values Msun ~ 2 • 10^" kg, apioneer ~ 8 • 10" ^° 

ms~^ and Gjy « 7 • 10~^^ m^ s~^ kg, we find 

k = 2'k/L K, i ■ 10" ^^m"^. This provides an estimate 
for L: 



L w 2 • 10 



(9) 



One light year being equal to 9, 461 • 10^^ m, we find that 
L is more or less one fifth of a light year. This justifies, a 
posteriori, why the Pioneer anomaly is so weak while the 
usual Newtonian theory is so efficient at the scale of the 
solar system. 

These results are in agreement with the constraints de- 
rived so far, which are: 

(1) L differs from 0. 

(2) L is quite larger than the size of the solar system. 
Actually, although a more detailed treatment is out of 

the scope of the present paper, it is not too difficult to 
convince oneself that, in the approximation where we treat 
the sun as a sphere where p™"**'^'' is taken to be constant 



■^Up to the so-called Pioneer anomalies [8] that we discuss soon, 
and up to relativistic corrections which, anyhow, remain out of the 
scope of the present work, merely aiming at describing the motion 
of slowly moving objects relatively to "fixed stars". 



mate the perturbation resulting from (the modification of 
the usual Newton-Poisson equation by) new source terms 
including the vacuum density, as well as the back reac- 
tion of matter and vacuum resulting from the equivalence 
principle (as described at the level of equation ([5])). This 
refined model leads thus to the same prediction as before, 
which is, it predicts as main perturbation at the level of 
the solar system the existence of a constant acceleration 
directed towards the sun, of magnitude Gm ■ Msun -k^ 12. 

Coherent interferences in the far field regime. If we 
consider a distribution of matter of size smaller than or 
comparable to the size of the solar system (like a typical 
star for instance), the potentials generated by each ele- 
ment of the matter distribution will coherently interfere in 
the "far field regime" (|x — x'| >> L). A striking novelty 
of our approach, compared to the "standard" Newtonian 
approach is that, in this case, the gravitational field gen- 
erated by a mass M located at the position x', at the posi- 
tion X, is no longer the Newton-Coulomb potential "'^"^^ 



but it is the oscillating function 



-G.YjV/co.s(fc(|x-x'|) 



-x'l ' 

This 



leads to several counterintuitive properties, the first of 
them being that gravitational forces are not necessarily 
attractive. Indeed, the acceleration a^^'^^ generated by 
a mass located at position x' onto a test-particle located 
at position x obeys 



viND „ -g , -G'jvMcos(fc(|x- x'l) 
° ~ dr^ 47r|x-x'| ' 

-kGNMsin{k{\x - x'l) GNMcos{k{\x - x'|) 



(10) 



Therefore, in the far field regime (at distances larger than 
one light year), the dominating term of the potential 
created by any celestial object of size comparable to or 
smaller than the size of the solar system is, generically, 
the first one in the last part of equation pUj) : 



VIND 



-kG N M sin{kr) 



(11) 



The expression pT|) obviously differs from the expres- 
sion ([2]). Now, as we have shown in the previous section, 
L being smaller than one light year, it is significatively 
smaller than typical galactic sizes (even smaller than typ- 
ical sizes of galactic cores), which implies that the poten- 
tial generated by the stars of the galactic core will not 
always constructively interfere. As we shall now show, it 
is possible to find a good agreement between MOND and 
VIND models in this case (incoherent interferences in the 
far field regime) , provided we interpret the MOND theory 
in a statistical sense. 

Incoherent interferences in the far field regime. An- 
other counterintuitive feature of our model is that gravita- 
tional potentials generated by different sources of matter 
do not always add to each other in a simple way but have 
the ability to interfere with each other. This property is 
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obvious if we consider electro-magnetic potentials but is 
totally absent from the Newtonian paradigm. To illustrate 
the consequences of this property, let us assume that the 
core of a galaxy is densely populated so that the law of 
large numbers is valid. The acceleration generated by the 
core of the galaxy onto a test particle (or another star) 
located in the suburbian zone of the galaxy will then obey 
(where the symbol ...w... must be understood here as... is 
equal "in average", or is "typically" equal to...): 



VIND\ 



-C^/N^GNk < Mp, > 



o 



(12) 



where C is a constant of the order of 1 while < Mq > 
is the average mass of a star of the galactic core and A^0 
its number of stars. It must be so because the phases 
which characterize the potentials generated by the dif- 
ferent local clusters formed by the Nq stars of the core 
will add incoherently, so that the average "phase factor" 
< sin{kri) > obtained after summing over their 

incoherent contributions will be Gaussianly distributed 
around zero, with a variance a equal to C ■ ^/Nq, with 
C a number of the order of unitjH, as predicted by the law 
of large numbers. 

Consequently, oY^^^ , the (averaged over space) value 
of the acceleration towards the core of the galaxy, pre- 
dicted in the VIND approach, will be statistically dis- 
tributed around the value 0, with a spread of the order of 

a-kGN<MQ> y/NQGNk<MQ > 



as announced at the level 



--C- 

of equation ([T^ . 

We find in this case that the ratio between < 



-.VIND 



> 



and 



-.MOND 



obeys 



la^^^^l ^, I Gn < Mq > 
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Injecting in equation (|13|) the typical values Msun 
2 ■ lO^o kg, ao w 1,2 • 10-1° 
we find that 



ms-2 and fc w 3 • 10"^^ m^^ 



a 
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< Mq > 



2,5 



< Mp, > 



M,. 



(14) 

This demonstrates how, in a statistical sense, MOND 
phenomenology can be derived from the assumptions that 
we originally formulated about vacuum energy, which are 
the basis of our VIND model. 

That this ratio is not exactly equal to one is not amazing 
for two reasons: 

(1) MOND's acceleration merely provides an order of 
magnitude. Indeed, in the framework of MOND's phe- 
nomenology it is well-known that one free parameter must 



^Prom now on we shall make the natural choice C = a{sin{kr)) 



1/2 = a/ /g " sin^xdx/2-K) 



be adjusted, for each galaxy, in order to take account of 
the mass of the galaxy. 

(2) the ratio ^J^^^^^ which appears in equation is 

likely to increase due to various sources of decoherence and 
to decrease due to the presence of dark positive matter- 
energy. For instance, due to the spread of k one should 
expect that the effective MOND acceleration diminishes 
with the distance to the center of the galaxy, due to deco- 
herence at large distances. 

A fine study of the calibration factors which arc intro- 
duced galaxy by galaxy provide a possible test aimed at 
discriminating VIND and MOND theories. In particular 

the scaling in terms of ^/^j^^ is not the same in both 

approached. Other tests are discussed in the last section. 

Open questions, discussion and conclusion . 

Other possible experimental tests. We mentioned 
three possible tests which confirm or could confirm that 
our ideas are well-founded (Pioneer anomalies, flat galac- 
tic velocity distributions and fine study of the calibration 
parameters and velocity curves), for which no new brutto 
experimental data must be collected [11[S1[5]. Other tests 
could be realized in the futureas we discuss now. 

Although both theories predict that the direction of the 
acceleration undergone by a test mass (for instance a star 
in the suburbs of a galaxy) is directed towards the cen- 
tre of the galaxy, the VIND and MOND approaches differ 
by an essential feature: the intensity of this acceleration 
slowly varies over space in the MOND model while, in the 
VIND model, it undergoes a strong relative change (it can 
even change its sign) at a scale L (supposedly of the or- 
der of a fifth of light year). This implies another main 
difference between both theories: MOND theory is per se 
a deterministic theory while, due to the random distri- 
bution of constructive and destructive interference effects 
over space, VIND predictions are per se statistical. It is 
only in average that both theories lead to comparable pre- 
dictions. This suggests several approaches aimed at dis- 
criminating between both explanatory schemes. The first 
of them consists of checking the sign of the centripetal ac- 
celeration undergone by stars located at the edge of dense 
galaxies, but this is not an easy task for experimentalists 
because we are dealing here with tiny velocities and accel- 
erations. In the case of nearly circular orbits around the 
galactic core, the radial velocities are particularly small 
and vibrations along the radial axis predicted by the VIND 
model should be particularly slow. 

Now, even if the results confirm that the sign of the cen- 
tripetal acceleration is nearly always positive, this does 
not necessarily mean that the VIND approach must be 
rejected. The reason is that the distribution of the trajec- 
tories of stars, in the VIND approach, is not necessarily 
homogeneously spread over space. Indeed, it is reasonable 
to suppose that the zones where gravitation is repulsive 



''MOND's acceleration scales like (< Mq >)^''^, while VIND's 
acceleration scales like (< Mq >)^. 
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(roughly regions where < sinikr) > is negative) are un- 
stable because there the test particle undergoes a repul- 
sive gravitational interaction which can not be counterbal- 
anced by the centrifugal force. It is therefore likely that 
the regions where sin{kr) is negative are sparsely popu- 
lated (which is confirmed by the fact that a flat velocity 
distribution has been observed: when the gravitational 
force is repulsive the distance to the galactic core and the 
tangential velocity are a priori uncorrelated) . The latter 
prediction, in turn, is likely to be submitted to experi- 
mental observations because it would imply that inside a 
given galaxy populated regions of extent L would coexist 
with empty regions of comparable size. Actually, in order 
to derive the distribution of trajectories of stars of a given 
galaxy, in the VIND approach, one must ultimately de- 
scribe the full history of these trajectories which is a very 
difficult task. Nevertheless, it is very unlikely that this his- 
tory conspires to reproduce the deterministic predictions 
of MOND theory and therefore a measure of statistical 
fluctuations of the velocities of stars situated at the same 
distance of the core of a galaxy should suffice in principle 
to discriminate between MOND and VIND predictions. 

Remark that the VIND theory also predicts the appear- 
ance of repulsive barrers around a star, at distances sepa- 
rated by integer multiples of L, which could maybe be re- 
vealed by a direct study of accretion disks of stars present 
in media where a lot of interstellar gas is present. This 
would be the firt experimental evidence of negative (re- 
pulsive) gravitational interaction. 

Another situation of interest occurs when the law of 
large numbers is not valid, for instance when a small num- 
ber of stars or a small number of galactic clusters are 
present. This is interesting because it is known that in 
this regime the MOND approach does not fit well with 
observations. Let us consider for instance that three stars 
are aligned, and separated by a distance larger than L, 
and let us treat the star situated in the middle as a test 
particle. In the MOND approach it will undergo an ac- 
celeration ^GAr(Mi/r^3 — M2/rl^) ■ aq. In the VIND ap- 
proach however it will undergo an acceleration equal to 
I + G^M..m(fcri3) | Qbviously these predic- 

tions have few in common. 

Conclusion. Before we end this last section, it is worth 
noting that our model indirectly provides an estimate of 
the cosmological constant A = • A:^ (67r)^ • 10""'^^ 
4 • 10~^^s~^, to compare with the lower value of 10~'^^s~^, 
which has been derived on the grounds of rather recent 
astronomical observations [3] (aimed at measuring indirect 
manifestations of the existence of a non null cosmological 
constant predicted in the framework of general relativity), 
and with the larger value of 10®''s~^ which derives from 
a cut-off at the Planck scale [B]. Despite of the fact that 
it is not clear today, to anybody, what is the meaning 
of the cosmological constanIO and how and why it should 

^Averaging the acceleration expresssed in equation l)10|l over a 
spatial period, wc sec that gravitation vanishes. Rather ironically, 



be related to quantum vacuum fluctuations and quantum 
vacuum energy, we were pleased to elaborate on the idea 
that a negative and constant "etheric" distribution of mass 
is omnipresent in the cosmos. As we have shown, the 
recognition of this fact has implications at the level of 
astronomy which cannot be ignored, and it is not necessary 
to resort to sophisticated theories (like general relativity, 
exotic particles and so on) in order to figure out what some 
of these implications could be. 

The most amazing lesson that we can draw from our 
analysis is that (in the case where our theoretical building 
would appear to resist to facts) it allows us to infer "the 
weight of empty space" . Indeed, it predicts for the vac- 
uum density of mass —2px = —k'^c^ /AttGn a value close to 
— 0, 9-10"'^ kg/m"^, more or less MINUS one gram per cubic 
meter. ..Once we opened the Pandora box, many questions 
spontaneously arise, e.g., "is k constant throughout time 
?", "to which extent does its existence influence the cos- 
mologic evolution (expansion, galactogcncsis, distribution 
of the cosmic background and so on) ?", "can its value 
be infered from the standard particle model ?"...As we 
see, black clouds in the blue sky of astrophysics possibly 
generate rain and rainbows. 

* * * 

As an outsider to the field who essentially built his cul- 
ture of the topic through the reading of magazines, wiki 
sites (e.g. the site of John Baez dedicated to vacuum's 
energy (http://math.ucr.edu/home/baez/vacuum.html') , 
the wiki sites respectively dedicated to the cosmological 
constant and the MOND theory) and so on, the author 
wishes to apologize in advance in the case where the bib- 
liography would be incomplete. 
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